Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 23 November 2011 



(MN JATeH style file v2.2) 



Ring galaxies from off-centre collisions 



o 

(N 

> 
O 

(N 
(N 

o 

u 

6 

(N 
> 
oo 
(N 
O 
O 



M. Mapelli^ & L. Mayer^ 

^ INAF-Osservatorio astronomico di Padova, Vicolo dell'Osservatorio 5, 1-35122, Padova, Italy, michela.mapelli 
^ Institute for Theoretical Physics, University of Ziirich, Winterthurerstrasse 190, CH-8057, Ziirich, Switzerland 



inaf . it 



23 November 2011 



ABSTRACT 

We investigate the formation of RE galaxies (i.e. of coUisional ring galaxies with an 
empty ring), with N-body/SPH simulations. The simulations employ a recipe for star 
formation (SF) and feedback that has been shown to be crucial to produce realistic 
galaxies in a cosmological context. We show that RE galaxies can form via off-centre 
collisions (i.e. with a non-zero impact parameter), even for small inclination angles. 
The ring can be either a complete ring or an arc, depending on the initial conditions 
(especially on the impact parameter). In our simulations, the nucleus of the target 
galaxy is displaced from the dynamical centre of the galaxy and is buried within 
the ring, as a consequence of the off-centre collision. We find that the nucleus is not 
vertically displaced from the plane of the ring. We study the kinematics of the ring, 
finding agreement with the predictions by the analytic theory. The SF history of the 
simulated galaxies indicates that the interaction enhances the SF rate. We compare 
the results of our simulations with the observations of Arp 147, that is the prototype 
of RE galaxies. 
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1 INTRODUCTION 

Ring galaxies are characterized by a bright ring of gas and 
stars, whose diameter may be very large (< 100 kpc, Ghosh 
& Mapelli 2008) and where star formation (SF) is generally 
very intense (Higdon 1995, 1996; Marston & Appleton 1995; 
Mayya et al. 2005; Romano, Mayya & Vorobyov 2008). A 
large fraction of ring galaxies (« 60 per cent. Few & Madore 
1986) is thought to form via (almost) head-on collisions with 
massive intruder galaxies (at least one tenth of the target 
mass). Because of the gravitational perturbation induced by 
the bullet galaxy, a density wave propagates through the disc 
of the target galaxy, generating an expanding ring of gas and 
stars (e.g., Lynds & Toomre 1976; Theys & Spiegel 1976; 
Toomre 1978; Appleton & James 1990; Appleton &i Struck- 
Marcell 1987a, 1987b; Struck-Marcell & Lotan 1990; Hern- 
quist & Weil 1993; Mihos & Hernquist 1994; Appleton & 
Struck-Marcell 1996; Gerber, Lamb & Balsara 1996; Struck 
1997; Horellou & Combes 2001; Mapelh et al. 2008a, 2008b). 
Therefore, ring galaxies are unique laboratories where the ef- 
fects of galaxy interactions can be studied, from a plethora 
of different points of view. First, the kinematics of the ring 
provides information about the dynamics of the interaction. 
The simple geometry of ring galaxies, combined with the 
essential kinematic data, allows us to reconstruct the main 
features of the interaction and the time elapsed since it oc- 



curred. Furthermore, ring galaxies are often characterized by 
high SF rate (SFR, up to « 20 M© yr"^), suggesting that 
the density wave associated with the propagating ring trig- 
gers the formation of stars. When the density wave reaches 
the external region of the target galaxy, SF may involve 
'fresh' gas, almost unpolluted by previous episodes of SF. 
This likely explains why the only three ring galaxies (the 
Cartwheel, AMI 159-530 and Arp 284) for which metallicity 
measurements have been published (Fosbury & Hawarden 
1977; Weilbacher, Due & Fritze-v. Alvensleben 2003; Smith, 
Struck & Nowak 2005) have metallicity Z < 0.3 Zq (when 
recalibrated on the basis of the electron temperature or of 
the P- method, see, e.g., Pilyugin, Vflchez & Thuan 2010). 

One particular subclass of ring galaxies, named RE by 
Theys & Spiegel (1976), is characterized by the fact that the 
ring is empty at its interior, apparently lacking the nucleus. 
This class includes objects that are among the closest and 
the most studied ring galaxies (e.g., Arp 146, Arp 147, VII 
Zw 466, I Zw 28, NGC 985, AM 0058-220, AM 2145-543). 
Arp 147 is generally considered the prototype of this class. 
Arp 147 has a high SFR (~ 4 - 12 M© yr"\ Romano et al. 
2008; Rappaport et al. 2010; Fogarty et al. 2011), and a reb 
atively young ring: ~ 50 Myr elapsed after the collision, as 
estimated on the basis of the ring expansion velocity (113 ±4 
km s~^, Fogarty et al. 2011) and diameter (~ 15 kpc, de 
Vaucouleurs et al. 1991). Arp 147 has a total stellar mass 
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similar to the Milky Way (6 - 10 x 10^° Mq, Fogarty et al. 
2011). A fraction ~ 0.04-0.34 of the total stellar light is con- 
tributed by a young (1 — 70 Myr) stellar population (Fogarty 
et al. 2011). Recently, Rappaport et al. (2010) reported the 
detection of 9 ultra-luminous X-ray sources (ULXs, i.e. non- 
nuclear point-like X-ray sources with luminosity Lx ^ 10^^ 
erg s~^) in the ring of Arp 147. 

Lynds & Toomre (1976) proposed that an off-centre col- 
lision can displace the nucleus and make it appear embed- 
ded in the ring. This idea is supported by observations (e.g., 
Thompson & Theys 1978; Schultz et al. 1991), that show 
that the suspected remnant nucleus can be distinguished 
from the other knots of the ring for its lack of SF. 

Gerber, Lamb & Balsara (1992) confirm, by running N- 
body/smoothed particle hydrodynamics (SPH) simulations, 
that RE galaxies form via off-centre collisions. They propose 
that the ring of RE galaxies is incomplete and may appear a 
complete ring because of the viewing angle. They find that 
the remnant nucleus is not only off-set from the centre, but 
also out of the plane of the ring, and appears to be buried 
within the ring only because of projection effects. On the 
other hand, Fogarty et al. (2011) find no evidence for the 
offset of the nucleus from the plane of the disc, based on 
kinematics. 

In this paper, we run N-body/SPH simulations of off- 
centre galaxy collisions, to study the formation of RE galax- 
ies. The simulations are described in Section 2. Section 3 re- 
ports the results and Section 4 summarizes the main conclu- 
sions. In our simulations, we adopt initial parameters close 
to the properties of Arp 147. However, we are not interested 
in producing a 'perfect' model of Arp 147, but in under- 
standing the formation of RE galaxies, from a more general 
point of view. 



2 SIMULATIONS 

The initial conditions for both the target and the intruder 
galax:y are generated by using an upgraded version of the 
code described in Widrow, Pym & Dubinski (2008; see also 
Kuijkcn & Dubinski 1995 and Widrow & Dubinski 2005). 
The code generates self- consistent disc-bulge-halo galaxy 
models, derived from explicit distribution functions for each 
component, that are very close to equilibrium. In particu- 
lar, the halo is modelled as a Navarro, Frenk & White (1996, 
NFW) profile. The disc (when present) is exponential. The 
code also includes the possibility of adding a Hernquist bulge 
(Hernquist 1993). We generate bulgeless models for the tar- 
get galaxy and models with bulge for the bullet galaxy. 

Masses and characteristics lengths of the simulated 
galaxies are listed in Tables 1 and 2. We report the results of 
eight runs. Runs A, B, C and D will be considered our fidu- 
cial models, whereas runs E, F, G and H were performed for 
comparison with Gerber et al. (1992). The only differences 
among runs A, B, C and D are the initial mass of gas and 
stars of the target galaxy. The total initial mass in baryons 
(gas and stars) is the same (8.48 x 10^'^ M0) for all the sim- 
ulations, but the ratio between gas and star mass (/gas) is 
different: /gas = 0.5, 0.2, 0.09 and 0.05, in run A, B, C and 
D, respectively. These values of /gas cover the typical range 
observed for Milky Way-like galaxies (see, e.g., Haynes et al. 
1999; Geha et al. 2006). 



Runs E, F, G and H have the same star and gas content 
as run D. They differ from run D for the impact parameter 
(runs E and G) and/or for the target-to-companion mass 
ratio (runs F, G and H). Run H differs from the other runs 
also for the halo-to-baryon mass ratio (~ 2.5 instead of ~ 8, 
for comparison with Gerber et al. 1992) and for the zero 
inclination angle. 

We set the mass of the companion to be approximately 
half of (equal to) the mass of the target in runs A, B, C, D 
and E (F, G and H).^ 

We initially set the centres of the target and the bul- 
let galaxy at a distance of three virial radii of the target. 
For runs A, B, C, D, F and H, the impact parameter is 
6 = 8 kpc. In runs E and G b is = 12 kpc. In all runs, 
the initial relative velocity between the centres of mass is 
Vrei = 900 km s"^ (see MapeUi et al. 2008a). The incUna^ 
tion angle (with respect to the symmetry axis of the target) 
is Q = 7°, with the exception of run H (where a = for 
comparison with Gerber et al. 1992). The inclination an- 
gle and the relative velocity are slightly different from those 
obtained for Arp 147 by the most recent measurements^. 
In the Appendix A, we show a simulation performed with 
Q = 33°, in agreement with the results by Fogarty et al. 
(2011). However, the morphology of the ring formed in this 
simulation looks very different from that of Arp 147, as too 
much warping is induced by such a large inclination angle. 
We will investigate the role of inclination angle and warping 
in a forthcoming paper. Furthermore, we stress again that 
the aim of this paper is to investigate the formation of RE 
galaxies, rather than producing a perfect model of Arp 147. 

The mass resolution of the simulations is 10^ M© for 
stars and gas, and 5 x 10* M© for dark matter (DM). The 
softening length is 0.1 kpc. Wc simulate the evolution of the 
models with the N-body/SPH tree code gasoline (Wadsley, 
Quinn & Stadel 2004). Radiative cooling, SF and supernova 



^ Romano et al. (2008) report a mass ratio Mtarg/Mint =0.57 
between the ring galaxy and the intruder in the Arp 147 sys- 
tem. They derive the two masses from the K magnitude and 
from the colours (using the recipes by Bell et al. 2003). Ac- 
cording to their estimates the stellar mass of the ring galaxy 
Arp 147 is A/, = 2.11 X 10^°Mq. On the other hand, Fogarty 
et al. (2011) estimate a stellar mass of the ring galaxy Arp 147 
M« = 6.1 — 10.1 X 10^" Mq, by combining optical measurements 
and galaxy population models. Since Bell et al. (2003) models are 
known to be less accurate in the case of late-type galaxies than 
in the case of early-type galaxies, it is reasonable to adopt the 
mass estimate by Romano et al. (2008) for the intruder galaxy 
and that by Fogarty et al. (2011) for the target. Therefore, the 
mass ratio between target and intruder is Mtarg/A^int = 1.6 — 2.7 
(assuming that the dynamical mass scales approximately as the 
stellar mass). Consequently, we adopt Mtarg/Mint = 2 in runs 
A, B, C, D and E. We also maJce some check runs (F, G and 
H), where the mass of the target is the same as the mass of the 
intruder, for comparison with Gerber et al. (1992). 
^ Fogarty et al. (2011) constrain the line-of-sight maximum dis- 
tance between the two galaxies to be ( = 7.8 kpc and measure a 
projected distance of d = 12.8 kpc. They also estimate the empty 
ring to be inclined by j = 25° with respect to the plane of the 
sky. Prom simple geometrical arguments, this leads to a minimum 
inclination angle for the interaction of ~ 33° , larger than our fidu- 
cial model (in our model, a = 7° and Urei = 900 km imply 
/ = 30 kpc, d=15 kpc and i = 20°). 
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Table 1. Initial conditions. 





Target 


Bullet 


M,'' [10^° Mq] 


see Table 2 


1.91 


Gas Mass'' [10^° Mq] 


see Table 2 




Halo scale length'^ [kpc] 


6.0 


6.0 


Disc scale length [kpc] 


3.7 


3.0 


Disc scale height [kpc] 


0.37 


0.30 


Bulge scale length [kpc] 




0.6 



M* is the total stellar mass of the galaxy. The stars of the 
target are distributed according to an exponential disc. The 
stars of the bullet are distributed according to an exponential 
disc (75 per cent of the total stellar mass) and an Hernquist 

bulge (25 per cent).'' The gas of the target is distributed 
according to an exponential disc, with the same parameters 
(scale length and height) as the stellar disc. We name halo 
scale length the NFW scale radius Rb = R200/C, where -R200 is 
the virial radius of the halo (NFW 1996) and c the 
concentration (here we assume c = 12 for both galaxies). 




Figure 1. Projected mass density of stars and gas in run C. 
The simulated target galaxy has been rotated by 20° along the 
x-a,xis in the direction of the y-axis and by 20° counterclockwise, 
to match the observations of Arp 147. The scale is logarithmic, 
ranging from 11.15 to 7.04 X 10* Mq pc~^. From top to bottom 
and from left to right: 10, 30, 50 and 70 Myr after the interaction. 



(SN) blastwave feedback are enabled, as described in Stinson 
et al. (2006, see also Katz 1992). The adopted parameters 
for SF and feedback (see Section 3.4) are the same as used in 
recent cosmological simulations capable of forming realistic 
galaxies in a wide range of masses (e.g., Governato et al. 
2010; Guedes et al. 2011). 



3 RESULTS 

3.1 Morphology of the ring in runs A, B, C and D 

We first consider the morphological evolution of the ring in 
our four fiducial runs (A, B, C and D). Fig. 1 shows the 
time evolution of run C. The ring starts developing 20 — 30 
Myr after the collision. About 50 Myr after the collision the 
diameter of the ring is 15 — 19 kpc (similar to the observed 
diameter of Arp 147). We notice that the target galaxy de- 
velops a complete ring. 

Another interesting feature of Fig. 1 is the structure of 
the companion galaxy. We simulate the bullet as an early- 
type spiral galaxy (with a purely stellar disc and without 
gas) . The initial conditions of the bullet were chosen to allow 
for the formation of a bar, sufficiently strong to produce a 
ring instability in the disc of the bullet (e.g., Buta & Combes 
1996; Tiret & Combes 2007). We evolved the companion in 
isolation till the formation of the resonant ring and we adopt 
this configuration as initial condition. This choice was made 
to test whether the companion of Arp 147 is a resonant ring 
galaxy. The resonant ring of the bullet survives relatively 
unperturbed for the first ^ 50 Myr after the interaction, 
but it is more and more perturbed in the later stages of the 
evolution. This result cannot be used to constrain the prop- 
erties of the companion of Arp 147, since we consider only 
one possible model for the companion. However, this may in- 
dicate either that the resonant ring of Arp 147 companion is 
more stable than our model against the perturbation (e.g., 
because of a more massive halo), or that we observe the 
Arp 147 system ^ 50 Myr after the collision. Both cases 
are allowed by the data (e.g., Fogarty et al. 2011). A fur- 
ther possibility is that even the ring of the companion has 
a collisional origin (that is, it was formed in the same inter- 
action as the ring of the target). This hypothesis was pro- 
posed, e.g., for II Hz 4 (Lynds & Toomre 1976). Although we 
cannot exclude this scenario for Arp 147, there are various 
hints against it. First, the primary is an empty ring, pro- 
duced by a collision with a large impact parameter, whereas 
the secondary is a symmetric ring, with a centred nucleus 
(although the system is observed with a large inclination), 
consistent with a nearly axisymmetric collision. These two 
morphologies are unlikely to form as a consequence of the 
same interaction. Second, the ring of the secondary is very 
smooth and regular, as observed in resonant ring galaxies 
(e.g., NGC 3081, Buta & Purcell 1998). 

Figs. 2 and 3 show the projected density of gas and 
stars of the target (seen face-on), respectively, in runs A, B, 
C and D. The gas has a clumpy structure and is confined 
in a thinner ring than the stellar component. Both gas and 
stars form complete rings. In the simulations with a higher 
fraction of gas (runs A and B), the distribution of stars is 
clumpier, because the young stellar population is significant, 
and tracks the distribution of gas. Instead, in the runs with a 
lower fraction of gas (runs C and D), the old stars dominate 
the mass distribution. 

Fig. 4 shows the surface density profile of gas (top) and 
stars (bottom) in runs A, B, C and D (which differ only for 
the gas to star mass fraction) . The four runs behave in a very 
similar way, with a strong density peak in correspondence 
of the ring (radius r 6 — 10 kpc, at t = 50 Myr after the 
collision) . 

Fig. 5 shows the projected density of gas (top) and 
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Table 2. Differences in the initial conditions between runs. 









Run A 


Run B 


Run C 


Run D 


Run E 


Run F 


Run G 


Run H 


Target 


DM Mass [10" 


Mq] 


7.0 


7.0 


7.0 


7.0 


7.0 


7.0 


7.0 


2.1 




M. [IQiO Mq] 




5.66 


7.06 


7.77 


8.13 


8.13 


8.13 


8.13 


8.13 




Gas Mass [10^° 


Mq] 


2.82 


1.42 


0.71 


0.35 


0.35 


0.35 


0.35 


0.35 


Bullet 


DM Mass [10" 


Mq] 


3.7 


3.7 


3.7 


3.7 


3.7 


7.4 


7.4 


2.2 


Impact Parameter b [kpc] 




8 


8 


8 


8 


12 


8 


12 


8 


Inclination a \°] 




7 


7 


7 


7 


7 


7 


7 






In this Table, we show only the parameters that have been changed in different runs. 





Figure 2. Mass density of gas of the target galaxy, projected in 
the X — y plane, at t = 50 Myr after the collision. From left to 
right and from top to bottom: run A, B, C and D. The scale is 
logarithmic, ranging from 2.22 to 1.40 X 10^ Mq pc'^. 



Figure 3. Mass density of stars of the target galaxy, projected 
in the x — y plane, at t = 50 Myr after the collision. From left to 
right and from top to bottom: run A, B, C and D. The scale is 
logarithmic, ranging from 2.22 to 2.22 X 10* Mq pc'^. 



stars (bottom) when the target galaxy is seen edge-on. Only 
runs A and D are shown, but runs B and C are intermediate 
cases between the plotted ones. We find that, in all the con- 
sidered runs, the nucleus of the target galaxy is still inside 
the disc even in the ring phase: the collision has off-set the 
nucleus from the dynamical centre and buried it in the ring, 
but the nucleus is not vertically displaced out of the ring. 
Actually, the entire ring, including the nucleus, is vertically 
displaced (by 4 kpc) with respect to the initial plane of 
the target disc. 

3.2 Morphology of the ring in runs E, F, G and 
H: comparison with Gerber et al. (1992) 

In our runs A, B, C and D, both the stellar and the gas ring 
are complete. This is not the case of Gerber et al. (1992) 
runs: they conclude from their simulations that Arp 147 ring 
is incomplete and appears complete only because of projec- 



tion effects. This discrepancy is likely due to differences in 
both the initial conditions and the resolution between our 
simulations and that by Gerber et al. (1992). In particular, 
the main differences between our runs A, B, C and D and 
those by Gerber et al. (1992) are the following, (i) Gerber 
et al. (1992) assume a mass ratio Mtarg/Mint =1 between 
target and intruder, respectively, while in the runs shown 
above we assume Mtarg/A/int =2. (ii) There is no DM halo 
in the galaxies simulated by Gerber et al. (1992), but just a 
spherical stellar halo (modelled as a King profile), whereas 
we simulate a NFW DM halo. Furthermore, the halo to disc 
mass ratio in Gerber et al. (1992) is only ~ 2.5. (iii) The 
mass and spatial resolution of Gerber et al. (1992) are a fac- 
tor of > 10 worse than those of our simulations, (iv) In our 
runs A, B, C and D, the impact parameter is fo = 8 kpc, cor- 
responding to a distance = 4 — 5 kpc between the centres 
of mass of the two galaxies during the collision. In the paper 
by Gerber et al. (1992), D = 7 kpc, which means that the 
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Figure 4. Top panel: surface density profile of gas of the target 
galaxy in run A (dotted line, red on the web), run B (solid line, 
blue on the web), run C (dashed line, green on the web) and 
run D (black dot-dashed line), at t = 50 Myr after the collision. 
Bottom panel: surface density profile of stars of the target galaxy, 
at t = 50 Myr after the collision. Lines are the same as in the top 
panel. 



collision is much more off-centre with respect to our simula- 
tions, (v) The simulations by Gerber et al. (1992) have no 
inclination angle. 

To check the importance of these differences, we per- 
form runs E, F, G and H, exploring various parameters. In 
particular, runs E and G have a larger impact parameter 
(b = 12 kpc), runs F, G and H have Mtarg/Mint =1. Run H 
was set to reproduce most of the features of Gerber et al. 
(1992) simulations (i.e., a DM to baryon fraction ~ 2.5, 
Mtarg/Mint =1, no inclination angle, D = 7 kpc). 

Fig. 6 shows the density map of gas (face-on) in sim- 
ulations E, F, G and H, at f = 50 Myr after the interac- 
tion. The ring is incomplete in these four runs, especially in 
run G (with Mtarg/Mnt =1 and & = 12 kpc, bottom left- 
hand panel). In runs E and G (top and bottom left-hand 
panel) the partial ring is also strongly deformed, as a con- 
sequence of the large impact parameter, combined with the 
non-zero inclination angle (which induces some degree of 
warping) . 

The ring in run H, although incomplete, is much more 
circular and less deformed than the other rings: this is due 
to the absence of warping in runs with zero inclination an- 
gle. The fact that the halo-to-disc mass ratio is smaller in 
run H than in the other runs has no significant effect on the 
propagation of the ring, at least in the initial stages. This 
is a consequence of the fact that the formation of the ring 
occurs in the inner regions of the target, that are baryon- 
dominated (we remind that the baryon mass in run H is 
the same as in the other runs). Finally, the ring in run H is 



Figure 5. Top panels: mass density of gas of the target galaxy, 
projected in the y — z plane, at t = 50 Myr after the collision. 
From left to right: run A and D. Bottom panels: mass density of 
stars of the target galaxy, projected in the y — z plane, at t = 50 
Myr after the collision. From left to right: run A and D. The scale 
is logarithmic, ranging from 7.04 X 10~^ to 7.04 X 10* Mq pc~^. 




Figure 6. Mass density of gas of the target galaxy, projected in 
the X — y plane, at t = 50 Myr after the collision. From left to 
right and from top to bottom: run E, F, G and H. The scale is 
logarithmic, ranging from 2.22 to 1.40 X 10^ Mq pc'^. 

smaller (by a factor of ~ 1.5) than the ring in run F, at the 
same time. This is likely due to the fact that the mass of the 
intruder is smaller (by a factor of ~ 3) in run H with respect 
to run F"^. In fact, we know from the analytic model that the 
radial velocity perturbation, in the impulse approximation, 

^ The total mass of the target is also smaller by a factor of ~ 3 
in run H with respect to run F, but the mass of the target in the 
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Figure 7. Mass density of stars of the target galaxy, projected 
in the y — z plane, at t = 50 Myr after the collision. From left to 
right and from top to bottom: run E, F, G and H. The scale is 
logarithmic, ranging from 7.04 X IQ-^ to 7.04 X 10* Mq pc'^. 



scales with the mass of the intruder (see, e.g., equation 2 of 
Struck-Marcell & Lotan 1990). 

In summary, a larger impact parameter is crucial to pro- 
duce incomplete rather than complete rings, in combination 
with other parameters (e.g., the target-to-intruder mass ra- 
tio). Therefore, we suggest that RE galaxies form partly as 
complete and partly as incomplete rings, depending on the 
impact parameter, on the involved galaxy masses, potentials 
and scale-lengths. 

From the analysis of runs A, B, C and D, we concluded 
that the nucleus of the target galaxy is not vertically dis- 
placed out of the ring (the entire ring, including the nucleus, 
is vertically displaced hy ^ 4 kpc with respect to the ini- 
tial plane of the target disc. Fig. 5). Instead, Gerber et al. 
(1992) propose that the nucleus is displaced out of the plane 
defined by the ring by ~ 3 kpc and it appears buried in the 
ring because of projection effects. 

Fig. 7 shows the edge-on projected density of stars in 
runs E, F, G and H. The edge-on view of the ring in these 
runs is more perturbed than in our fiducial runs (A, B, C 
and D). In particular, the disc is much thicker in those sim- 
ulations where Mtarg/Mint =1 (runs F, G and H). However, 
we cannot conclude that the nucleus is completely displaced 
with respect to the rest of the ring. We cannot make more 
accurate comparisons, as Gerber et al. (1992) do not provide 
more information about the displacement they find. We ex- 
pect that the possible differences between our simulations 
and those by Gerber et al. (1992) are connected with the 
different resolution or with the stability of the initial condi- 
tions. 

As we discussed in Section 2, we decided to simulate 
the target as a bulgeless galaxy. However, we also run a 
check case in which the target has a bulge and all the other 
properties are the same as in run D (see Appendix B). In 
particular, we want to check whether the nucleus of the tar- 
get can be more easily displaced from the plane of the ring 
after the collision, if the target has a bulge. The results in- 
dicate that the presence of a bulge (instead of a nucleus) in 



inner regions is substantially unchanged, as the baryon mass is 
unchanged. 



the target does not change the features of the interaction 
(e.g., the completeness of the ring, the position of the nu- 
cleus, the kinematics and SFR, etc.). In particular, there is 
no displacement of the nucleus from the plane of the ring 
after the collision. 

3.3 Kinematics 

The simulations also provide information about the kine- 
matics of the interaction. Fig. 8 shows the most relevant 
features of the velocity field of the simulated RE galaxy at 
t — 50 Myr after the collision (run A). The left-hand panel 
shows the tangential velocity «tan, i.e. the two-dimensional 
tangential velocity with respect to the disc of the target. 
Therefore, Vmn represents the rotational component of the 
target disc. The target galaxy rotates almost regularly for 
radii r 10 kpc, whereas the rotation is strongly perturbed 
for r < 10 kpc, because of the expanding ring. 

The central panel of Fig. 8 shows the radial velocity 
frad with respect to the disc of the target, and measures the 
expansion of the ring. Vra.d has various interesting features. 
First, the plot of Wrad highlights the high-velocity expanding 
ring at r ~ 10 kpc. The ring has a maximum of the expansion 
velocity for the negative a;— axis (which corresponds to the 
southern part of Arp 147). 

Second, the matter in the outer parts of the disc (r > 10 
kpc) in the region of the negative y— axis (which corresponds 
to the eastern part of Arp 147) has negative values of v^^d 
and is falling toward the centre. This means that the densest 
part of the ring (which corresponds to the negative j/— axis, 
see Figs. 2 and 3) is composed of both inner particles that 
are directed outwards and outer particles that are directed 
inwards. This result agrees with the predictions of the ana- 
lytic caustic theory (see, e.g., Struck-Marcell & Lotan 1990; 
Appleton & Struck-Marcell 1996; Struck 2010). 

Third, the matter in the outer parts of the disc (r > 10 
kpc) in the region of the positive t/— axis has Vrad ~ 0: it 
appears still unperturbed by the collision. Actually, this part 
of the target galaxy is the farthest from the position of the 
impact with the intruder. The right-hand panel of Fig. 8 
shows the ratio between |i7rad| and |utan| and highlights the 
strong expansion of the ring, that dominates locally over the 
rotation. 

Table 3 shows the average values of |utan| and l^radl 
in the inner part of the target galaxy (r < 10 kpc). These 
values are very similar for all the runs, as the orbits and 
the involved masses are almost the same. In particular, 
(|frad|) ~ 140 ± 80 km s"^ and (jwtani) ~ 200 ± 110 km 
s~^ (with the marginal exception of runs G and H). The 
radial velocity {|vrad|) is comparable to the value of the ex- 
pansion velocity derived by Fogarty et al. (2011) for Arp 147 
(Vcxp ~ 113±4 km s~^), whereas (lutanl) is quite higher than 
the observed rotational velocity (Kot = 47 ± 8 km s~^). 
However, the value of (Ivtanl) obtained from the simulations 
cannot be directly compared with the observations of Kot, 
because these two quantities are defined and derived in sub- 
stantially different ways. We also stress that the relative ve- 
locity between the simulated galaxies is quite different from 
the value obtained by Fogarty et al. (2011). Furthermore, 
the uncertainty associated with (|utan|) and with (|wrad|) is 
so large (Table 3) that these are consistent with the obser- 
vational results by Fogarty et al. (2011) for Arp 147. 
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Figure 8. Velocity field of gas of the target galaxy in run A, at t = 50 Myr after the collision. Left-hand panel; tangential velocity of 
gas (associated with the rotation of the disc) with respect to the disk of the target galaxy, projected in the x — y plane. Central panel: 
radial velocity of gas (associated with the expansion or infall of the ring) with respect to the disk of the target galaxy, projected in the 
X — y plane. Right-hand panel: ratio between the modulus of the radial velocity and the modulus of the tangential velocity, projected in 
the X — y plane. The scale is linear and color coding is indicated in the panels. 

peratures allowed by the adopted cooling function (a few 
thousand K). In the blastwave feedback model of Stinson et 
al. (2006) , the gas heated by the SN energy is not allowed to 
cool for a timescale of the order of 10 Myr, the exact value 
being self-consistently determined by the blastwave solution 
of the McKee & Ostriker (1977) model, that assumes a blast- 
wave produced by the collective effect of many type II SNae. 
For type I SNae, the energy is instead radiated away on the 
cooling time (since they explode on longer timescales they 
will explode at significantly different time and would hardly 
produce a collective blastwave). 

To assess the effects of the SFR in our simulations, we 
will focus on runs A, B, C and D, that differ only for the 
gas-to-star fraction. The top panel of Fig. 9 shows the SFR 
as a function of time in the four considered simulations. In 
all the simulations, the SFR has approximately the same 
trend; it increases sensibly after the galaxy interaction. The 
main difference among the simulations is the normalization 
of the SFR, that spans from ~ 8 Mq yr~^ in run D to ~ 300 
Mq yr^^ in run A. Observations of Arp 147 (Rappaport et 
al. 2010; Fogarty et al. 2011) indicate a SFR- 4 - 12 Mq 
yr~^, consistent with runs C and D. 

The central and bottom panels of Fig. 9 show the time 
evolution of the total stellar mass of the target (M,) and 
of the mass fraction of young stars (My/M*; we define as 
young stars those stars that have an age fago < 100 Myr at 
t = 50 Myr after the interaction). The total stellar mass of 
the four simulations, at f = 50 Myr after the interaction, 
spans from 7 x 10^° to 8 x 10^° M©. Therefore, aU these 
simulations are consistent with the stellar mass of Arp 147 
(found to be M, = 6.1 - 10.1 x 10^" Mq by Fogarty et al. 
2011). 

However, there are large differences in the role of the 
young stellar population between different runs. In fact, 
My/M« ranges from 5 X 10"^ (in run D) to 0.17 (in run A) 
at t = 50 Myr. Since Fogarty et al. (2011) find that the mass 
fraction of young stars in Arp 147 is < 1 per cent, we confirm 



Table 3. Kinematics at t = 50 Myr after the interaction. 



(|fradl>[kms 1] (|i;tan|>[km s 



Run A 


134 ± 82 


186 ± 117 


Run B 


135 ± 79 


187 ± 115 


Run C 


144 ± 82 


207 ± 115 


Run D 


149 ± 86 


218 ±115 


Run E 


140 ± 104 


216 ±114 


Run F 


143 ± 95 


156 ± 128 


Run G 


174 ± 100 


217 ± 130 


Run H 


114 ±87 


193 ± 103 



3.4 Star formation 

In the simulations, the SF and feedback recipes of Stinson 
et al. (2006) are adopted. Three parameters characterize the 
SF and feedback recipe; (a) the SF threshold nsF, (b) the SF 
efficiency esp, and (c) the fraction of SN energy that couples 
to the interstellar medium (ISM) esN- SF occurs when cold 
(r < 3 X 10^ K), virialized gas reaches a threshold density 
nsF = 5 atoms cm~'^ and is part of a converging fiow. It 
proceeds at a rate 

dp, Pgas 1.5 /^N 

-TT = eSF ■- OC Pgas (1) 

at tdyn 

(i.e. locally enforcing a Schmidt law), where p, and pgas 
are the stellar and gas densities, and tdyn is the local dy- 
namical time. We choose esF ~ 0.1. The relatively high SF 
density threshold that we adopt is the same as that used in 
recent cosmological simulations that form realistic galaxies 
(Guedes et al. 2011) and it has been shown to have a key 
role in determining a realistic inhomogeneous structure of 
the ISM (see, e.g., Mayer 2011). Its value is determined by 
enforcing that at least one SPH kernel (32 particles) is con- 
tained within a spherical volume of diameter equal to the 
local Jeans length for gas at density nsp at the lowest tem- 
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Figure 9. Top panel: SFR as a function of time in run A (open 
squares, red on the web), run B (filled circles, blue on the web), 
run C (filled triangles, green on the web) and run D (black 
crosses). Central panel: total mass in stars (M«) of the target 
galaxy. Symbols are the same as in top panel. Bottom panel: to- 
tal mass of young stars (My ) normalized to the total mass in stars 
of the target galaxy. Symbols are the same as in top panel. Time 
t = corresponds to the collision. 



that only runs C and D are consistent with the observations 
of Arp 147. 

The difference in My/M* among various runs is con- 
sistent with the difference expected from the Schmidt law. 
For example, there is a factor of ~ 34 between run A and 
D. From the Schmidt law, we expect that the SFR scales 
with Pgas, which means that there should be at least a fac- 
tor of ~ 23 difference between My in runs A and D. Since 
the mass of old stars in run D is a factor of ~ 1.4 higher 
than in run A, we expect a difference by a factor of ~ 32 in 
My/M, between run A and D. 



4 CONCLUSIONS 

In this paper, we investigate the formation of RE galaxies 
(i.e. of coUisional ring galaxies with an empty ring) with 
N-body/SPH simulations. We have shown that RE galaxies 
can form via off-centre collisions (i.e. with a non-zero impact 
parameter), even for small inclination angles (~ — 7°). 

In runs A, B, C and D (where fo = 8 kpc and 
Aftarg/Afint =2), the ring is always complete and it is not 
produced by a projection effect. In runs with a larger im- 
pact parameter (6 = 12 kpc) or with Mtarg/Mint =1, the 
ring looks incomplete, as proposed by Gerber et al. (1992). 
Therefore, we conclude that RE galaxies can have either 



complete or incomplete rings, depending on the parameters 
of the interaction. 

We find that the nucleus is displaced from the dynam- 
ical centre of the galaxy, as a consequence of the off-centre 
collision. However, the nucleus remains in the (perturbed) 
plane of the ring. The entire ring (together with the nucleus) 
is slightly ( ^ 4 kpc) vertically displaced with respect to the 
previous plane of the disc of the target galaxy. This result 
is different from the findings by Gerber et al. (1992), who 
claim that the nucleus appears buried within the ring only 
as a consequence of projection effects. We expect that this 
discrepancy is related to the lower resolution of Gerber et 
al. (1992). 

The kinematics of the ring was studied in detail, con- 
firming the predictions by the analytic caustic theory (e.g., 
Struck-Marcell & Lotan 1990). We compare our simulations 
with the observations of Arp 147, finding similar expansion 
velocity in the ring. 

Finally, we study the SF history of the RE galaxy, 
finding that the encounter enhances the SFR of the tar- 
get galaxy. Runs C and D match the observational prop- 
erties of Arp 147 (Rappaport et al. 2010; Fogarty et al. 
2011), having average SFR~ 10 - 20 M0 yr~\ total star 
mass M, = 8 X 10^° Mq and mass fraction of young stars 
My/M, ~ (0.5 - 1.6) X 10"^ 

We stress that we adopt a recent implementation of the 
SFR and of SNae in our simulations, whose results are in 
good agreement with observations (see. e.g., Stinson et al. 
2006, 2009). However, we do not model the chemical evolu- 
tion of the gas, which may be very interesting in ring galax- 
ies, as pointed out by recent papers (e.g., Bournaud et al. 
2007; Mapelh et al. 2009, 2010; Michel-Dansac et al. 2010). 
Future N-body models of RE galaxies (and, in general, of 
coUisional ring galaxies) should account for multi-phase gas 
(to describe the evolution of molecular and atomic gas) and 
include metallicity evolution, to shed light on the SF history 
of ring galaxies. 
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is logarithmic, ranging from 11.15 to 7.04 X 10* Mq pc'^. 
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APPENDIX A: LARGE INCLINATION ANGLE 

In this appendix, we show a simulation where the inclina- 
tion angle is a = 33° (in agreement with the estimates for 
Arp 147 by Fogarty et al. 2011). The velocity components 
were rearranged so that the impact parameter is still 8 kpc. 
All the other initial conditions of this simulation are the 
same as in run D. From Fig. Al, it is evident that the ring 
of the simulated galaxy is strongly warped, because of the 
large inclination angle. We will investigate in detail the ef- 
fects of warping in a forthcoming paper. 



APPENDIX B: TARGET GALAXY WITH 
BULGE 

In this appendix, we discuss the results of a simulation in 
which the target galaxy has a bulge (with a total mass of 
1.2 X 10^" Mq). The other properties of this simulation are 
the same as those of run D. Fig. Bl shows the projected 
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Figure Bl. Mass density of gas and stars of tlie target galaxy in 
the check run with bulge, at t = 50 Myr after the collision. Top 
panels: mass density of gas (left-hand panel) and of stars (right- 
hand panel) projected in the x — y plane (i.e., face-on). The scale 
is logarithmic, ranging from 2.22 to 1.40 X 10^ Mq pc-^ and from 
2.22 to 2.22 X lO" Mq pc'^ in the left-hand and in the right-hand 
panel, respectively. Bottom panels: the same as the top panels, 
but the target galaxy is projected in the y — z plane (i.e., edge-on). 
The scale is logarithmic, ranging from 7.04 X 10~^ to 7.04 X 10^ 
M0 pc-2. 

density of the target galaxy face-on (top) and edge-on (bot- 
tom) for both the gaseous (left-hand) and the stellar com- 
ponent (right-hand panel). No significant differences can be 
observed with respect to run D in Figs. 2, 3 and 5. In the 
case of the target with bulge, the ring is still complete, 
and the nucleus/bulge is buried within the ring, but it is 
not vertically displaced. The SFR (6 M© yr"^ at t = 50 
Myr after the collision) is also very similar to run D (8 Mq 
yr ~ ). Furthermore, the kinematics of the ring cannot be 
distinguished from run D, as in the simulation with bulge 
i^rad = 149 ± 89 km s'^ and vt^.n = 201 ± 111 km s'^ (we 
remind that Urad = 149 ± 86 km s~^ and t^tan ~ 218 ± 115 
km s~^ in run D, see Table 3). Therefore, the assumption 
about the bulge of the target does not significantly affect 
our results, in the case of very late-type galaxies, such as 
Arp 147. 



